ITERATIVE PROCESS FOR TDMA RADIOMOBILE COMMUNICATIONS 



DESCRIPTION 

Technical field 

The purpose of this invention is an iterative 
process for TDMA (Time Division Multiple Access) 
radiomobile communications . 
5 It may be applied to the European GSM radiomobile 

system [1] and its development to GSM+ + or to the 
American D-AMPS system. 

State of prior art 

10 TDMA systems can be used to distribute several 

users communicating through the same radioelectric 
channel, in time [1,3]. Each user is separated from the 
others by periodically assigning a timeslot to him in 
which he can transmit, and another timeslot in which he 

15 can receive. Protective intervals are allowed at each 
end of each timeslot in order to guard against 
|| imperfections in time synchronization between 
transmitters and interference between consecutive 
timeslots resulting from the plurality of pathways. 

2 0 Figure 1 attached thus shows timeslots TA 1 , TA 2 , 

. . . , TAy allocated to N users with protective intervals 
IG. 

The radioelectric signal transmitted during a 
timeslot is obtained by transposition of the equivalent 
25 base band signal into frequency. The base band signal 
is the result of filtering the data stream to be 
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transmitted through a "transmission" filter. Normally, 
the transmitted data stream is composed of two 
sequences of data symbols separated in time by a 
sequence of symbols known to the receiver, called the 
reference sequence. The two sequences of data symbols 
may possibly originate from coding and interlacing of 
information to be transmitted. 

The reference sequence frequently has good time 
correlation properties, in order to enable a precise 
estimate of the channel at the receiver. For example, 
the reference sequences used in the GSM system, called 
CAZAC (Constant Amplitude - Zero Autocorrelation) 
sequences [1, 4, 5] are sequences with components taken 
from the bipolar alphabet {-1,1} and possessing a null 
circular autocorrelation function everywhere except at 
the origin. 

The radiomobile channel used during a 
communication between a transmitter and a receiver is 
usually of the multipath type with fast fading called 
RAYLEIGH fading. The existence of several paths is due 
to the fact that the radioelectric wave is propagated 
along several paths between the transmission location 
and the reception location. 

The received signal is then the sum of several 
more or less delayed replicas, more or less modified in 
phase and in amplitude. In order to obtain a reliable 
reproduction of the information carried by these 
replicas, the receiver applies filtering matched to the 
transmission filter and to the channel and it combines 
energy contributions from all transmitted signal 
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replicas in an optimum manner. The output signal from 
the matched filter is sampled at the rate of the 
symbols and it is whitened using a discrete filter 
called a whitener. 

Samples at the output from the whitener filter 
supply a filtered and noisy version of the transmitted 
data stream. The discrete filter associated with this 
filtered version, also called the discrete channel, has 
a finite pulse response that varies from one sample to 
the other. It characterizes the production of the 
multipath channel during the corresponding reception 
timeslot, in an indirect manner. 

Samples at the output from the whitener filter 
corresponding to the reference sequence are used to 
estimate the discrete channel [1, 4, 5]. This estimate 
of the discrete channel is used to equalize the 
remaining samples, so that the two transmitted data 
sequences can be detected if there are any. The 
equalizer usually used is known as the VITERBI 
equalizer with flexible decisions [2]. As its name 
suggests, this equalizer uses an application of the 
Soft-Output Viterbi Algorithm (SOVA) [2] to produce soft 
decisions on all transmitted data symbols. These soft 
outputs are possibly de-interlaced and decoded to 
detect the information transmitted. 

Figure 2 attached shows these 

transmission/reception operations. The transmitter E 
comprises a data source 10, an encoder /interlacer 12 
outputting symbols a k , and a modulator 14. The 
multipath channel is symbolized by a block 20. The 
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receiver R comprises a demodulator 3 0 (matched 
f ilter/whitener filter) outputting samples R k , a 
discrete channel estimator 32, an equalizer 34, a de- 
interlacer/decoder 36 and finally an addressee 38. 

The discrete channel seen at the output from the 
whitener filter may vary significantly from one 
timeslot to the next. This variation is due mainly to 
the change in propagation conditions between the 
transmitter and the receiver and to the frequency 
stability at the receiver. 

Propagation conditions have a direct influence on 
the observed multipath channel. They change due to a 
modification in the environment or a displacement of 
the transmitter and/or the receiver. They create a 
time variation in the discrete channel, both between 
successive timeslots and within the same timeslot. 

The variation of the discrete channel between two 
successive timeslots allocated to the same user is 
particularly large when the time interval between these 
timeslots is large. This variation is accentuated, 
even at a single timeslot, due to an increase in the 
carrier frequency or the speed of the transmitter 
and/or the receiver. 

In practice, the variation of the discrete channel 
between two timeslots is sufficiently large to prevent 
any modulated estimate of the channel. The estimate of 
the channel in a timeslot must then be based solely on 
the samples of the corresponding reference sequence. 

The use of this reference sequence provides an 
unbiased estimate of the discrete channel at the mid- 
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point of the timeslot. Like the GSM system, an 
invariable channel is usually assumed within a received 
timeslot. In this precise case, the estimate obtained 
at the mid-point of the frame can be used without 
degradation for equalization of the rest of the frame. 
But high speed movements of some terminals and the 
sustained demand for higher rate services operating at 
increasingly high radioelectric frequencies, make this 
assumption less and less justified. The discrete 
channel may be affected by significant variations 
between the beginning and the end of a given timeslot. 
The difference between this real discrete channel and 
its estimate becomes increasingly large as the distance 
from the reference sequence increases. This may cause 
a large and irreversible deterioration in the reception 
quality and/or performances of the TDMA system. 

One solution [13] has already been proposed to 
solve this problem of fast variation of the channel in 
the case of radioelectric channels related to satellite 
radiocommunications . However, this solution only deals 
with the case of one channel with a single path and 
assumes that reference or data symbols have a common 
transmitted energy. Therefore, it cannot be applied to 
multipath channels encountered in land radiomobile 
communications . 

Furthermore, like the EDGE extension in the GSM 
system, high rates may be obtained with unchanged 
frequency band by the use of modulations using an 
increasing number of states. The system then becomes 
very sensitive to the channel estimating quality, and 
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6 

severe deterioration may occur even at very low 
displacement speeds if the power and/or length of the 
reference sequence are not increased. 

The purpose of this invention is to overcome these 
5 disadvantages. 



Disclosure of the invention 




The purpose of this invention is to improve the 
performances of TDMA systems by reducing the ratio 

10 between the received power and the multiple access 
interference for a given reception quality, or in other 
terms, to improve the reception quality for a constant 
reception power. This improvement in quality makes it 
possible to increase the capacity and coverage of TDMA 

15 systems such as GSM system, and particularly its EDGE 
extension to offer high rate services. In particular, 
this improvement is obtained due to an optimization of 
the channel estimate at the receiver. This 
optimization is a means of overcoming any deterioration 

2 0 in the performances generated by a fast variation in 
the channel between the beginning and the end of the 
same timeslot. It can also significantly reduce 

deterioration caused by an increase in the number of 
modulation states, without increasing the length or 

2 5 power of the reference sequence. 

Another purpose of this invention is to reduce the 
cost price of terminals by making them less sensitive 
to the poor frequency stability of an inexpensive local 
oscillator (for the same reception quality) , for 
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example such that the received signal can be transposed 
into base band. 

The invention can also reduce the number and/or 
power of reference symbols transmitted during each 
timeslot, for the same reception quality. This purpose 
is achieved by taking account of data symbols in the 
discrete channel estimate that are usually more 
numerous than the reference symbols. This purpose is 
also achieved by the optimum usage of samples 
corresponding to the reference symbols and/or data 
symbols of an arbitrary number of timeslots allocated 
to the user. This usage is achieved even in the case 
of large variations in the channel between two 
successive timeslots belonging to the same user. 
Finally, this purpose is achieved by the optimum usage 
of samples from the common signal channel, and possibly 
users on the down link. 

According to the invention, the reference symbols 
are distributed into groups throughout the duration of 
the allocated timeslot, which guarantees better 
tracking of the channel variation and therefore 
increased robustness with regard to high displacement 
speeds and poor frequency stability of the local 
oscillator. Reference sequences may be chosen 

independently of each other, but they must always have 
good autocorrelation properties. They may be derived 
from short CAZAC type sequences. 

With the invention, the energy transmitted per 
symbol may vary from one symbol to another, without 
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causing difficulties for the channel estimate which is 
always done optimally. 

According to the invention, a timeslot is 
processed by block every time that the samples 
corresponding to this timeslot and possibly to other 
timeslots (originating from the same user, from other 
users or from the common signal channel) are available 
to make the channel estimate. 

Like a conventional TDMA receiver (such as the GSM 
system receiver) , the first step is always a coarse 
estimate of the discrete channel using samples 
depending only on reference symbols only in the 
timeslot to be processed. This coarse estimate 

characterizes the phase and amplitude variation of all 
coefficients (improperly called "paths") in the 
discrete channel selected for carrying out the 
equalization, for each reference or data symbol in the 
block to be processed, and symbol by symbol. 

One of the characteristics of the invention is the 
use of a mesh to represent the multipath channel. This 
characteristic reduces the complexity of the receiver. 
If we consider the example of sequences of data symbols 
comprising (ji) symbols without coding, with a two-state 
phase modulation (MDP2), each transmitted bit may be 
equal to either of two values. Since these two values 
are equally probable, this means that there are 2 n 
possible sequences of data symbols. Therefore the 
complexity of a traditional calculation using partial 
probabilities for all possible sequences increases 
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exponentially (2 n ) with the number of symbols per 
sequence . 

The use of a mesh according to the invention to 
represent the channel considerably reduces the 
complexity of the calculations, since the calculations 
will be sub-divided into sections, each section 
requiring a certain calculation type, always the same 
for each section. Therefore, the complexity increases 
with the number of sections, which increases in 
proportion to the number n. Therefore the complexity 
is linear, rather than exponential. 

The labels of the branches of mesh sections change 
with the variation of the discrete channel with time in 
each of the timeslots forming the block. Furthermore, 
an algorithm other than the Soft-Output Viterbi 
Algorithm (SOVA) [2] is used for a more precise 
characterization of the probabilities of the branches 
of the mesh dependent on the samples received. If the 
channel estimator is used optimally, these conditional 
probabilities may be determined using a BAHL or BCJR 
(abbreviation for Bahl Coke Jelinek Raviv) algorithm 
[6]. Other simpler algorithms [10], usually derived from 
linearization of this algorithm, may also be used, 
however performances will be slightly reduced. They 
can considerably reduce the complexity of the channel 
estimator . 

In the case of a conventional receiver, the only 
weighted outputs provided by the SOVA algorithm are 
used directly after de-interlacing to decode the 
transmitted information. In the case of the receiver 



SP 16348. C RS 



10 



according to the invention, the conditional 
probabilities indirectly characterizing the data 
symbols are used in addition to the reference symbols 
to provide a better quality estimate of the discrete 
channel . 

The improved estimate of the discrete channel 
obtained at the end of a given iteration may be used to 
improve the quality and reliability of conditional 
probabilities of branches of the mesh characterizing 
the discrete channel. These improved conditional 

probabilities may then be used jointly with the 
reference symbols to provide an additional improvement 
to the estimate of the discrete channel. 

In theory, an iterative estimate will not give an 
optimum estimate of the discrete channel until after an 
infinite number of iterations has been carried out. 
However in practice, a few iterations are sufficient to 
obtain performances that are almost as good as an 
optimum estimate. 

Once the discrete channel estimate has been made, 
the corresponding mesh may be used either by the BAHL 
optimal algorithm (or by any simplification of it) , or 
by the SOVA algorithm to provide improved weighted 
outputs to the decoder. 

The optimum nature of the receiver according to 
the invention is related to the quality of the discrete 
channel estimate seen during a given timeslot. This 
optimum nature is based on the use of an iterative 
algorithm called the SAGE (Space-Alternating 
Generalized Expectation-Maximization Algorithm) [7] to 
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find the most probable channel set up according to 
samples received from the block to be processed. The 
SAGE algorithm is an extension of the EM (Expectation- 
Maximization) [8,9] algorithm, that eliminates the 
5 problem of coupling between coefficients of the 
discrete channel while they are being estimated. 

The estimate of the discrete channel is also based 
on the decomposition of each discrete channel path to 
be estimated according to a KARHUNEN-LOEVE [3] expansion 
10 algorithm. This decomposition provides firstly a soft 
characterization of the variations with time of each of 
the discrete paths, and secondly is naturally 
integrated into the SAGE algorithm. It can also 
process cases in which the energy transmitted per 
15 symbol varies from one symbol to another. 

Finally, the discrete channel estimate is based on 
the use of the BAHL algorithm which provides the 
probabilities of the mesh branches characterizing the 
discrete channel to each iteration of the SAGE 
20 algorithm, according to received samples and the 
channel estimate supplied in the previous iteration. 

Therefore more precisely, the purpose of the 
invention is a radiomobile communication process of the 
time division multiple access type, in which timeslots 
25 are allocated to several users for transmission and 
reception of radioelectric signals, and in which: 

- on transmission, data symbols and reference 
symbols are created in each allocated timeslot, 
a radioelectric signal containing these data 
3 0 symbols and reference symbols is transmitted, 



SP 16348. C RS 



12 



this signal being sent along a multi-path 
radiomobile channel , 

- on reception, matched filtering is applied to 
the received signal, the filtered signal is 
sampled at the rate of the symbols, an estimate 
of the radiomobile channel is made, the samples 
are processed using this channel estimate and 
data symbols specific to the allocated timeslot 
are restored, 

s process being characterized in that: 

- on transmission, several reference sequences 
possessing appropriate time correlation 
properties are built up for each allocated 
timeslot, and these sequences are distributed in 
each timeslot, 

- on reception: 

i) in order to estimate the radiomobile 
channel, this channel is represented by a 
mesh comprising branches symbolizing 
transitions between two successive states of 
the channel and characterizing possible 
sequences at the channel output, 

ii) the probabilities of these branches are 
calculated using an iterative algorithm in 
which an estimate of a representation 

{g 1(D) of L paths of this channel is 

calculated using a finite number (D+l) of 
iterations using the maximum a posteriori 
probability criterion, the iteration rank 
(d+l) (3d+l) being used to obtain a re- 



16348. C RS 



13 



estimate {G' (rf+1) } L 1 starting from the estimate 

{G l(d) JJ" ' obtained in the previous iteration 
(3d) , and in each iteration (3d+l) , the 
vectors {G ,(d+1) } 1=0, 1, L-1 are 

calculated in L steps from vectors {G m(d+1) }' m ^ 

and |G m(d) ) L 1 calculated in the previous 

step, the estimate {c 1(D) Jj* ' of L channel 
paths finally being obtained from the 
estimate of its representation {g I(D) } . 

Brief description of the drawings 

• Figure 1, already described, shows the 
distribution of timeslots allocated to several 
users ; 

• Figure 2, already described, is a diagram 
illustrating transmission and reception using 
the TDMA technique ; 

• Figure 3 shows variations of the bit error rate 
(BER) as a function of the normalized Doppler 
spreading (B D T S ) ; 

• Figure 4 shows variations of the bit error rate 
(BER) as a function of the Es/No ratio for MDP4 
and MDP8 modulations, in a first configuration 
and for a discrete channel with two paths; 
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• Figure 5 shows variations of the bit error rate 

(BER) as a function of the E5/N0 ratio for the 
MDP4 and MDP8 modulations in a second 
configuration with a discrete channel with two 
5 paths ; 

• Figure 6 shows variations of the bit error rate 

(BER) as a function of the E5/N0 ratio for an 
MDP4 modulation using the second configuration 
and a discrete channel with four paths; 
10 • Figure 7 shows variations of the bit error rate 

(BER) as a function of the Es/N ratio for an 
MDP4 modulation using a discrete channel with 
one, two and four paths; 

• Figure 8 is a view of the discrete channel seen 
15 at the output from the whitener filter, in the 

form of an offset register; 

• Figure 9 illustrates the estimate of the filter 
coefficients of the discrete channel by 
correlation of received samples with a CAZAC 

20 sequence and its circular offsets; 

• Figure 10 shows the distribution of the largest 
eigenvalues of the covariance matrix for three 
normalized Doppler spreadings; 

• Figure 11 illustrates a representation of the 
25 mesh of the normalized discrete channel seen at 

the output from the whitener filter for a 
channel with two paths and a modulation with two 
states and a length of the timeslot equal to 
four; 
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• Figure 12 shows the chronology of processing 
done on reception for a received timeslot; 

• Figure 13 shows the chronology of processing 
done in an elementary iteration of the channel 

5 estimating algorithm. 

Description of particular embodiments 

The perf ormances of the process according to the 
invention have been compared with the performances of a 

10 variant of the conventional process used for the GMS 
system [1]. A conventional receiver in the GSM system 
starts by estimating the discrete channel at the mid- 
point of the timeslot to be processed by correlation of 
received samples that depend solely on the reference 

15 sequence with the corresponding CAZAC sequence. It 
then equalizes the other received samples using the 
SOVA algorithm in order to decode the transmitted 
information symbols. The SOVA sub-optimal algorithm 
was replaced by the BAHL optimal algorithm used in the 

20 invention, in order to equitably evaluate the intrinsic 
improvement in performances provided by the invention 
compared with a conventional GSM receiver. 

The comparison is made using the variation in the 
unmodified bit error rate (BER) (without considering 

25 any error correction or detection coding) , for 
reference symbols grouped at the mid-point of the 
timeslot to be processed (like the GSM system) , and 
distributed in small groups throughout the duration of 
this timeslot (according to the invention) in order to 
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give better tracking of the variation of the channel 
with time. 

Figure 3 thus shows the variation of the BER as a 
function of the normalized Doppler spreading B D T S , which 
5 is the product of the Doppler spreading B D and the 
symbol duration T s that characterizes the variation 
intensity of the channel due to displacement of the 
transmitter and/or receiver. Doppler spreading is 
proportional to the speed of the mobile terminal and to 

10 the carrier frequency. The frequency instability of 
the local oscillator, which is ignored here, introduces 
additional performance deteriorations . Like the GSM 
system, the behavior of the BER was evaluated for 
timeslots composed of N=148 symbols each. The MDP4 

15 modulation was selected (modulation with 4-state phase 
displacement) . The discrete channel considered is 
composed of two paths (L=2) with the same average 
power. The Es/M ratio between the average energy 
received per symbol and the noise level is 10 dB. 

2 0 Three configurations were considered: 

■ A central reference sequence consisting of 2 6 
bipolar symbols taken at random among the 8 
sequences used in the GSM system. These 
sequences are cyclic extensions to the 8 CAZAC 

25 sequences with length 16 [1, 11]. 

■ a central reference sequence 

-1, +1, +1, +1, -1 , of 5 bipolar symbols 
obtained by cyclic extension (except for a sign 
and an offset) of the unique bipolar CAZAC 
30 sequence with length 4, namely -1, +1, +1, +1 . 
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■ Five groups of reference sequences distributed 
across the entire timeslot to be processed, 
identical to the previous central reference 
sequence of the five bipolar symbols (total of 
25 reference symbols). 
Figure 3 shows representative curves corresponding 
to the following cases: 

41: conventional receiver, first configuration, 
42: conventional receiver, second configuration 
43: receiver according to the invention, second 

configuration, 
44: conventional receiver, third configuration, 
45: receiver according to the invention, third 
configuration . 

These results show that when sequences grouped at 
the mid-point of each timeslot (first and second 
configurations) are used, it is impossible to go 
outside normalized Doppler spreading B D T s =l/2500 without 
significantly deteriorating performances . Conversely, 
the use of groups of small distributed sequences (third 
configuration) makes it possible to achieve 
standardized Doppler spreading close to B D T S =1/100 
without any perceptible deterioration. This 
significant improvement in the robustness of the 
receiver with regard to the Doppler effect, and 
possibly the frequency instability of the local 
oscillator, is due to the distribution of reference 
sequences, which makes it possible to more precisely 
monitor any variation of the channel with time. 
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Table 1 contains figures characterizing the 
improved performances provided by the invention with 
regard to the improved conventional algorithm used in 
the GSM system. In this table, C.R. means 

"conventional receiver" and I.R. means receiver 



according to 


the invention . 


The Es/No ratio between 


the average 


energy received per symbol and 


the noise 


level is 10 


dB . The improvement due to the 


invention 


is negligible in the case of 


a long central 


reference 


sequence (first configuration) 






Normalized 


First 


Second 


Third 


DOPPLER 


configuration 


configuration 


configuration 


spreading 


C.R. I.R. 


C.R. I.R. 


C.R. I.R. 


(B D T S ) 








1/2500 


0.0099 0.00855 0 


.01584 0.0089 


0.0105 0.0077 


1/100 


X X 


X x 


0.0143 0.0083 




Table 


1 





In this case, the reference sequence is 
sufficiently long to guarantee a good channel 
estimating quality, without the use of unknown data 
sequences. However, in the case of a single short 
reference sequence (second configuration) , the receiver 
according to the invention can give better performances 
than a conventional receiver, but also the performances 
are comparable to performances obtained with the long 
reference sequence of 2 6 reference symbols. 
Consequently, the receiver according to the invention 
is made insensitive to even a large reduction in the 



SP 16348. C RS 



size of the reference sequence. This result is 

explained by the fact that the channel estimate makes 
use of symbols in data sequences that cooperate with 
reference symbols. In the case of five groups of short 
reference sequences (third configuration) , the receiver 
according to the invention is capable of keeping the 
bit error rate almost unchanged by varying the 
normalized Doppler spreading between 1/2500 and 1/100. 

Figures 4, 5 and 6 show the bit error rate as a 
function of the E s /No ratio between the average energy 
received per symbol and the noise level for normalized 
Doppler spreading B D T s =l/2500 which is not very severe. 

Like the GSM system, the behavior of the BER was 
evaluated for timeslots composed of N=148 symbols each. 
However, the selected modulations are MDP4 and MDP8 
modulations. The discrete channel considered is 

composed of one (L=l) or two (L=2) or four (L=4) paths 
with the same average power. Two configurations were 
considered: 

■ A central reference sequence of 2 6 bipolar 
symbols taken from the 8 sequences used in the 
GSM system. These sequences are cyclic 
extensions of the 8 CAZAC sequences with length 
16 [1, 11]. 

■ A central reference sequence -1, +1, +1, +1, -1 
for two paths and -1, +1, +1, +1/ -1, +1/ +1 
for four paths with 4+L-l bipolar symbols 
obtained by cyclic extension of the unique 
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(except for a sign and an offset) bipolar CAZAC 
sequence with length 4, namely -1, +1, +1, +1 

For figure 4, the selected modulations are MDP4 
and MDP8 modulations using the first configuration with 
two paths (L=2). The curves are applicable to the 
following cases: 

- 51: MDP8, conventional receiver, 

- 52: MDP8, receiver according to the invention, 

- 53: MDP8, receiver with perfectly known 

discrete channel, 

- 54: MDP4, conventional receiver, 

- 55: MDP4, receiver according to the invention, 

- 56: MDP4, receiver with perfectly known 

discrete channel. 



In figure 5, the selected modulations are the MDP4 
and MDP8 modulations using the second configuration and 
a discrete channel with two paths (L=2) . The curves 
are applicable to the following cases: 

- 61: MDP8 , conventional receiver, 

- 62: MDP8 , receiver according to the invention, 

- 63: MDP8, receiver with perfectly known 

discrete channel, 

- 64: MDP4, conventional receiver, 

- 65: MDP4, receiver according to the invention, 

- 66: MDP4, receiver with perfectly known 

discrete channel . 
In figure 6, the selected modulation is the MDP4 
modulation with the second configuration and a discrete 
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channel with four paths (L=4) . The curves are 

applicable to the following cases: 

- 71: conventional receiver, 

- 72: receiver according to the invention, 

- 73: receiver with perfectly known discrete 

channel . 

The curves in figures 4 and 5 show that, unlike a 
receiver according to the invention, the performances 
of a conventional receiver in the GSM system are very 
sensitive to the length of the reference sequence used 
and deteriorate significantly between the first 
configuration and the second configuration. The curves 
in figures 5 and 6 also show that performance 
deteriorations with a conventional receiver increase 
significantly with the number of paths although 
conversely these deteriorations remain almost unchanged 
and are low for a receiver according to the invention. 
Therefore, the receiver according to the invention has 
excellent robustness with respect to any reduction in 
the length of the reference sequence and any increase 
in the number of channel paths. 

Table 2 contains figures characterizing the 
improvement in performances provided by the receiver 
according to the invention compared with the improved 
GSM system algorithm, for the MDP4 modulation. The B D T S 
ratio is equal to 1/2500, and the bit error rate is 
equal to 10" 2 and N=148. 
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Gain in dB 


Number of 


Long 


Short 


Receiver according to the 


pa. ens v ij / 


■v- o d ~y~ c^y~\ r** a 


reference 


invention with short 




sequence 


sequence 


sequence , conventional 
receiver with long sequence 


1 


0.03 


1.03 


0.03 


2 


0.48 


1.83 


0.48 


4 


1.42 


3 .04 


1.43 






Table 


2 



This table shows that a receiver according to the 
invention always has better performances than an 
5 improved conventional receiver, for a given length 
(first or second configuration) of the reference 
sequence. These results also show that the improvement 
provided by the invention is particularly significant 
when the length of the reference sequence is short or 

10 when the number of paths in the channel is increased. 
These results also show that even with a short 
reference sequence (second configuration) , a receiver 
according to the invention always has better 
performances than a conventional receiver with a long 

15 reference sequence (first configuration). 

Figure 7 shows the bit error rate as a function of 
the Es/No ratio for severe normalized Doppler spreading 
B D T S =1/100 and timeslots each composed of N=100 symbols. 
2 0 The MDP4 modulation is always selected. The discrete 
channel considered is composed of one (L=l) or two 
(L=2) or four (L=4) paths with the same average power. 
A single configuration with three distributed reference 
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sequences is considered. Each sequence is composed of 
4+L-l bipolar symbols and is obtained by cyclic 
extension (expect for a sign and an offset) of the 
single bipolar CAZAC sequence with length 4, namely - 

5 1, +1, +1, +1 . 

The value of Doppler spreading considered is 
sufficiently severe so that a conventional receiver 
cannot follow the variations of each path of the 
discrete channel that may be affected by phase 
10 rotations equal to ±Jt between the beginning and end of 
each timeslot. 

The curves in figure 7 are applicable to the 
following cases: 

- 81: L=l, receiver with perfectly known discrete 
15 channel, 

- 82: L=l, receiver according to the invention 

- 83: L=2, receiver with perfectly known discrete 

channel 

- 84: L=2 , receiver according to the invention 

20 - 85: L=4, receiver with perfectly known discrete 

channel 

- 86: L=4, receiver according to the invention. 
The curves in figure 7 show that the receiver 

according to the invention has very little 
25 deterioration of performances compared with an ideal 
receiver that knows the discrete channel perfectly. 
This deterioration is of the order of 0.4 dB for one 
path (L=l) , 0.5 dB for two paths (L=2) and 1 dB for 
four paths (L=4) . 
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Now that we have described the performances of the 
process according to the invention, the reception phase 
and particularly the received signal processing will be 
described in more detail. 
5 A receiver according to the invention determines 

the production of the most probable discrete channel 
during a given timeslot, depending on the samples 
received from the block to be processed. This 
determination is made using a Maximum A Posteriori 

10 (MAP) probability criterion. This channel estimate 
requires knowledge of a priori probabilities of symbols 
in transmitted data sequences. 

In order to simplify the description, the 
following mathematical discussions and the 

15 corresponding figures are always applicable to the case 
in which each block to be processed for the channel 
estimate is composed solely of received samples related 
to a single timeslot. The case in which symbols 
(reference or other symbols) from other timeslots (from 

20 the same user or from other users, from the common 
signaling channel, etc.) are used, is dealt with in the 
same way. 

The transposition operator is denoted ( . ) ' . As 
shown in figure 2, every timeslot for a given user is 
25 composed of a series of N symbols a 0 , a 1# . . . , a N _ 1 that 
can be written in the form of a column vector 

a = ( a 0 , a 2 , . . . , a N _ 1 ) 
the N symbols possibly being phase modulated (MDP2, 
MDP4, MDP8, ...). The transmitted energy E k associated 
3 0 with each symbol a k is independent of values taken on by 
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this symbol but it may vary from one symbol in the 
timeslot to another. This variation is usually the 
consequence of the TDMA system allocating different 
energies to data sequences and to reference sequences . 
5 Also as shown in the same figure, samples at the output 
from the whitener filter corresponding to this timeslot 
form a column vector 

with N+L-l components, where L is the number of paths 

10 in the discrete channel. 

Figure 8a shows all means between the input to the 
transmitter modulator and the output from the receiver 
whitener filtering. Therefore this figure shows a view 
of the real radioelectric channel and the processing 

15 done in the receiver (demodulation, matched filtering 
and whitening filter) . This figure shows the modulator 
91, a block 92 symbolizing the multi-path channel, the 
demodulator 93 with its matched filter 94 and its 
whitener filter 95. 

20 Figure 8b shows a view of these means in the form 

of an offset register with (L-l) cells Y 1# Y 2 , Y L-1 . 

This register receives the current symbol a n and the 
cells contain the symbols with a delay of one symbol 
period namely a n _,, a n _ 2 , a n . (L-1) . The circuit shown 

25 also comprises (L) multipliers Z 0 , Z 1# . . . , Z L-l that 
contain (L) coefficients C^,C^,...C^ 1 respectively, an 
adder S connected to the multipliers, and a final adder 
receiving a signal w n representing a noise. This final 
adder outputs a signal r n . 
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Each received sample may be written in the 
following form: 

R k =X c L a k-i +w k ' k=0 ' 1 — N + L-2, 
1=0 

where is the 1 th coefficient of the discrete channel 
filter at the k th sampling instant and W k is a white 
Gaussian additive noise. To simplify the presentation 
and make concise equations, the coefficients for the 
discrete channel satisfying k < 1 or k > N+l are used 
in the expression of R k although they are equal to zero. 
For the same reason, the transmitted symbols a K 
satisfying k < 0 or k > N-l appear in this expression 
of although they are equal to zero. 

The discrete channel coefficients are almost 
decorrelated between different discrete paths 
(independent of the variable 1) . Ignoring the effects 
of modulation and demodulation, this property is due to 
the fact that the multiple paths in the real 
radioelectric channel are independent. However, the 
coefficients of the same discrete path 1 are correlated 
with each other (dependent on the variable 7c) . Their 
discrete autocorrelation function defined by 

<t>L„=Ek(c , „)J 

depends only on the Doppler power spectrum for the 
corresponding path. 

Doppler power spectra usually encountered in 
radiomobile communications are either conventional (in 
external environments) or flat (in internal 
environments) . The autocorrelation function for the I th 
path with average power <j>o is given by 
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(j) 1 ^^ J 0 (7ikB D T s ) 

in the case of a conventional Doppler power spectrum 
and by 

^[=^ 0 sin (7ckB D T s )/7tkB D T s 

5 in the case of a flat Doppler power spectrum, where 
J 0 (.) is the first type of Bessel function with order 
0 [3, 12]. The average power §q varies from one path to 
another and characterizes the multipath intensity 
profile [3] of the discrete channel. 
10 The receiver according to the invention must use 

the maximum a posteriori (MAP) criterion to determine 
the most probable implementation of the vectors 

c 1 = (cj , c|_, ,cj. 2 ,... f c 1 ^., ) , 1 = 0.1, ...,L - 1 
characterizing the L discrete paths selected by the 

15 receiver for equalization. Consequently, it needs an 
appropriate representation of the discrete channel that 
it can use to become independent from the time 
correlation of coefficients for the same path. It also 
needs a representation that enables it to be 

2 0 independent of the variation (if any) of transmitted 
energy between one symbol in the timeslot and the next. 
This representation is obtained firstly by adding the 
vector of normalized transmitted symbols 

A = ( A 0 f A 1 , . . • , A N _ 1 ) 

25 with components A k =a k l^jE^ with a constant modulus 
equal to one. It is also obtained by rewriting the 
samples received during a timeslot in the following 
form: 
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R k =Z C k A k-i + W k ,k=0,l,...,N + L-l, 
1=0 

where C[ is the k th component of the associated 
normalized vector 

C 1 = (cJ^/E 0 ,cj +1 ^E7,cJ +2 ^/E 2 ,...,c n+m -7e n _! ) 

and directly characterizing the I th path in the discrete 
channel. Finally, it is acquired using a KARHUNEN- 
LOEVE decomposition [3] of each normalized vector C* in 

an orthogonal base {b*}^ taking account both of the 
variation of energy in transmitted symbols and 
characteristics of the Doppler power spectrum for the 
I th path [3]. With this decomposition, the normalized 
vector C (with correlated time components) of the 1 th 
discrete path can be described by an equivalent column 
vector G 1 = (g^GJ,...^ 1 ^) , in which the N components are 
well decorrelated. The vector C may be built up from 
its appropriate representation G' by the following 
transformation 

k=0 

Conversely, the components of this appropriate 
representation G 1 may be found within the same vector C 1 
by the inverse transformation 

G^ =(B ,k )* 1 C , ,k = 0,1,...,N-1, 
where (.)* is the conjugation operator. 

In accordance with the properties of the KARHUNEN- 
LOEVE decomposition [3], the vectors of the orthonormal 
base associated with the I th path are simply the 
eigenvectors of the co-variance matrix H* of the 
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normalized vector C' associated with this path. 
Furthermore, the associated eigenvalues ^A, j k=0 are 
simply the variances of the G^Gj^.^G^., components of 
the appropriate representation G' of the corresponding 
path. Therefore, the vectors {B lk } k J associated with 
each path, and the corresponding eigenvalues can be 
obtained from the matrix equations 

H l B*=X ft B\k = 0,l,...,N-l 

The (m, n) th entry in the co-variance matrix H 1 of 
the normalized vector associated with the 2 th path is 
given directly by 

HL = E[c m (C 1 . ) *]= VE m E n (t>Ln . 
where <t> m _ n denotes the discrete autocorrelation function 
of this path defined above. 

The receiver according to the invention must 
normally have L orthonormal bases associated with each 
of the L paths in the discrete channel in order to 
rebuild or decompose any intermediate or final estimate 
of the discrete channel. In practice, the shapes of 
the Doppler power spectra for the L channel paths are 
identical. Therefore the associated co-variance 

matrices {h 1 }^ are identical, except for a 
multiplication factor, and the corresponding 
orthonormal bases are also identical . A single 

orthonormal base in the receiver is then usually 
sufficient to estimate the channel satisfactorily. 

In order to generate the orthonormal base for the 
2 th path, the receiver must have the most precise 
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possible estimate of the co-variance matrix H' for this 
path. Consequently, it must have the most precise 
possible knowledge of the characteristics of the 
Doppler power spectrum of the radioelectric channel. 
This knowledge is regularly updated and must include an 
estimate of Doppler spreading B D , or an upper limit to 
it. It may possibly comprise an estimate of the 
average power of each path in the discrete channel and 
the shape of the Doppler power spectrum. Depending on 
this knowledge, the receiver may either calculate 
appropriate orthonormal bases in real time, or search 
in a bank of pre-calculated orthonormal bases for the 
bases that best match the characteristics of the real 
channel . In general , the receiver does not have 
precise knowledge of the shape of the Doppler power 
spectrum, or perfect and instantaneous knowledge of 
Doppler spreading. In this typical case, it may be 
assumed that the Doppler power spectrum is flat, and 
that the real Doppler spreading is equal to the upper 
limit of Doppler spreading. Furthermore, the receiver 
does not need to have a precise estimate of the average 
power of each of the discrete channel paths taken 
individually. Also in this case, a discrete channel 
with L paths with equal average powers may be assumed. 

The estimate of the c' vectors characterizing the 
L discrete channel paths is made iteratively using the 
MAP criterion through the normalized vectors C or their 
appropriate representation G' . This estimate may take 
optimal account both of the characteristics of the data 
symbols and the transmitted reference symbols 
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(statistics and energy) and the characteristics of the 
multi-path channel . 

The MAP estimate jf^ 1 }^ of one embodiment of the 

appropriate representation {g'JIo of the discrete 
channel is equal to 

fc'ju =argmaxp({3 , }|' =0 , |R) , 

iP J=o 

which maximizes the conditional a posteriori 
probability density p(fe ! ]u|R) • According to the 

invention, a solution as close as necessary to the 
exact solution can be obtained by iteration using the 
SAGE extension [7] of the EM algorithm [8, 9]. 

In general, the conditional a posteriori 

probability density pCfe'l^R) to be maximized has 
several global maxima, frequently leading to an 
ambiguity in the channel estimate according to the MAP 
criterion. This ambiguity may be cleared using 

reference symbols known to the receiver. However, this 
if often insufficient, since this conditional 
probability density also has local maxima that can be 
reached by the SAGE algorithm instead of the single 
global maximum. In order to avoid this problem, 
reference symbols have to be used to appropriately 

determine the initial conditions {C^}^ . 

The SAGE algorithm successively re-estimates each 
of the G' vectors by induction, keeping the previous 
estimates of the other vectors G™ unchanged, with m * 1 . 
In accordance with the properties of the SAGE 
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algorithm, these successive re-estimates guarantee 
monotonous growth of the conditional a posteriori 

probability density p({G ! }[l 0 1 |R) . 

Given the vector R of samples received during a 
5 timeslot to be processed, the receiver begins by 

calculating the initial conditions {g I(0) }J1o of vectors 

{(yjjlo starting from samples that only depend on 
reference sequences. Consequently, the reference 

sequences must be such that a fast and simple 

10 separation of the contributions of the L paths to each 
of these samples is guaranteed. In order to be able to 
make an initial estimate to get the channel estimating 
algorithm started, one or several reference sequences 
are necessary each having good autocorrelation 

15 properties. Like the GSM system, these sequences may 
be built starting from CAZAC sequences by extension or 
circular swapping. These sequences are composed of 
symbols with a constant modulus and with a zero 
circular autocorrelation function everywhere except at 

20 the origin. For the purposes of this invention, it is 
preferable to have CAZAC sequences with symbols 
belonging to modulation constellation used to transfer 
data symbols. For example, for an MDP2 modulation, the 
shortest CAZAC sequence alone (except for a sign and a 

25 circular swap), namely {-1, +1, +1, +1} only has four 
symbols. For an MDP4 modulation, using the {1, i, -1, 
i} constellation, the shortest CAZAC sequence alone 
(still with a multiplication factor i,-l or -i, and 
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except for a circular offset), namely {l,i} is composed 
of only 2 symbols. 

A reference sequence based on a CAZAC sequence 
with length M can only estimate up to M paths of the 
discrete channel. It is obtained by cyclic extension 
of L-l positions of this CAZAC sequence and contains a 
total of M+L-l symbols. Therefore the length M of the 
CAZAC sequence needs to be greater than or equal to the 
number of paths L of the discrete channel to be 
estimated. For the purposes of this invention, it is 
sometimes desirable to use short CAZAC sequences for a 
small number of paths L to be estimated (less than or 
equal to 4 for an MDP2 modulation and 2 for MDP4 , MDP8 
modulations, etc.), since they can significantly reduce 
the total number of reference symbols per timeslot, 
while guaranteeing more regular and more frequent 
sampling of the discrete channel, particularly in an 
environment with large Doppler spreading. 

Figure 9 illustrates the estimate of a discrete 
channel composed of two paths (L=2) in the case of an 
MDP2 modulation. The first line A gives the reference 
sequence, the second line Tl corresponds to the first 
path, the third line T2 corresponds to the second path 
and the fourth line S gives the composition of the 
received samples. The elementary reference sequence 
used {-1/ +1/ +1/ +1/ -1} composed of 5 symbols (M+L- 
1=5) is simply the circular extension of a single 
position in the CAZAC sequence {-1, +1, +1/ +1} 
composed of four symbols (M=4) . As shown in figure 5, 
the elementary reference sequence occupies positions 



SP 16348. C RS 



between n-4 and n in the timeslot to be processed. The 
rece ived samples are dependent solely on this sequence, 
and occupy positions between n-3 and n. They are given 
explicitly by: 

R„-3 = C n-3 A n-3 +C n-3 A n-4 + W n-3 =C - ^.3 + W„. 3 
Rn-2 = C°. 2 A.-2 + C'„. 2 A-3 + Wn-2 = C°. 2 + C' n . 2 + Wn-2 

R_ =C°.,A n _, +C' n . 2 A n _ 2 + W„, = C°, +C' n , +W n ., 
R n = C°A n +W n = -C° +C[ +W n 

If the estimate of the coefficients C£ and C' k in 
the normalized discrete channel is to be only slightly 
biased for a value of k between n-3 and n, it has to be 
assumed that Doppler spreading is sufficiently low to 
guarantee negligible channel variation between 
positions n-3 and n. In this case, an approximate 
estimate of these coefficients is obtained by 
correlating the samples, where k is between n-3 and 
n, with the CAZAC sequence itself and a circular offset 
of a position towards its left. The correlation with 
the CAZAC sequence {-1, +1, +1, +1} gives: 

CL =C' n . 2 =C' n , =C' n =(-R n _ 3 +R n _ 2 +R n _, +R n )/4 = 

(cU-HC , n . 2 +cL+c 1 n )/4 + w: 

as a first estimate of the samples n-3, n-2, n-1 and n 
in the second path, where is a term that takes 

account of noise and bias. The correlation with the 
sequence obtained after circular offset {+1, +1, +1, 
-1}, gives 

Cl> =C° n . 2 = C°, = C° a = (R n . 3 + R n _ 2 + R + R n )/ 4 = 
(C 3+ C°. 2+ C°. 1+ C:)/4 + W n ° 
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as the first estimate of samples at instants n-3, n-2 , 
n-1 and n of the first path, where W° is a term that 
plays the same role as . 

In the general case, the coarse estimate of the L 
5 channel coefficients at a given instant in the timeslot 
is obtained by correlation of samples that depend 
solely on the reference sequence with L offset versions 
of the corresponding CAZAC sequence. 

Each reference sequence built up by circular 

10 extension of a CAZAC sequence can be used to obtain a 
noisy sample of each of the discrete channel paths. 
The number of samples per timeslot necessary for 
reconstruction of the channel depends on the normalized 
Doppler spreading B D T S and the number N of symbols 

15 making up this timeslot. As the product B D TgN 

increases, the variation in the channel from one end of 
the timeslot to the other increases, and the number of 
elementary reference sequences needs to be greater. An 
elementary reference sequence is a sub-set of grouped 

20 symbols in the reference sequence associated with the 
timeslot obtained starting from a circular extension of 
a CAZAC sequence. 

In accordance with the properties of the KARHUNEN- 
LOEVE decomposition of the discrete channel, a higher 

2 5 value of the eigenvalue X 1 * for a given path 1 will lead 
to a higher average value of the contribution G^B 1 * of 
the eigenvector associated with the appropriate 
representation of this path. Figure 10 shows the 
distribution of the largest eigenvalues for a timeslot 

30 of N=100 symbols with the same energy normalized to 1 
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and normalized Doppler spreading B D T S equal to 1/50 
(curve 113), 1/100 (curve 112), and 1/500 (curve 111) 
respectively. This figure clearly shows that an 

appropriate representation of each of the paths of the 
discrete channel can be based on the few eigenvectors 
that have the largest eigenvalues, with no significant 
loss of performance. For low values of the product 
BpTgN, the eigenvector with the largest eigenvalue alone 
is sufficient to give a precise representation of each 
of the paths of the discrete channel. This eigenvector 
then varies very slowly with time, and therefore a 
reference sequence grouped at the mid-point of each 
timeslot can be used to rebuild a discrete channel that 
is practically constant over the entire timeslot. This 
is the case for the GSM system in which normalized 
Doppler spreading is assumed to be very weak and the 
channel is assumed to be constant over a given 
timeslot . 

For slightly higher normalized Doppler spreading, 
the two eigenvectors associated with the two largest 
eigenvalues are necessary for a precise reconstruction 
of the discrete channel. These two vectors vary 
sufficiently so that a reference sequence composed of 
two elementary sequences located on each side of the 
mid-point of the timeslot are necessary. 

In the general case, as the value of the B D T S N 
product increases, the number of eigenvectors that 
needs to be used in the appropriate representation of 
the discrete channel also needs to be increased. 
Furthermore, as this product increases, the variations 
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in these eigenvectors also increase, and therefore the 
number of elementary reference sequences to be used 
also increases. Furthermore, for a constant length of 
the reference sequence, as the number of elementary 
reference sequences increases, the length of each 
sequence needs to be reduced. For a given number of 
paths L, the length of an elementary sequence M+L-l can 
be reduced to 2L-1 at best, by choosing a CAZAC 
sequence with length L; this is why it is useful to 
have very short CAZAC sequences . According to the 
invention, the CAZAC sequence {-1, +1, +1, +1} with 
length M=4 is recommended for all phase modulations 
when the number of paths L is less than or equal to 4 
and the product R D TJ$ is low. Also according to the 
invention, the CAZAC sequence {l,i} with length M=2 is 
recommended for a number of paths less than or equal to 
2 and for fairly large products B D TJ<S. 

Let K be the number of elementary reference 
sequences per timeslot to be processed. Let 
a-A/ + i,d-A/ + 2,...,q, 1= 0,1,... ,L-1, be the discrete channel 
samples procured by the Jc th elementary reference 
sequence and S k = {k - M + 1 , k - M + 2 , . . . , k} be the 
set of indices of these samples. The initial 

conditions jp l(0) }!lo ' of the appropriate channel 
representation according to the invention may be 
obtained from these samples using the following 
equations : 

X S C nV B n k ) k _, N . 

of ,. w ;--L or •«» 

IN k m=0neS m ^ k 0 
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The weighting factor W k ' depends on the normal Nj. 
of the restriction of eigenvectors to include only 
samples corresponding directly to channel samples 
procured by the reference sequence. It also depends on 
the quality of the contribution of the eigenvector B lk 
measured by the X lk /N Q ratio. As the variance X lk of the 
coefficient in the representation of the discrete 

channel increases in comparison with the variance N 0 of 
noise samples, the corresponding weighting also 
increases and approaches 1. 

Still based on the vector R of samples received 
during a timeslot, the receiver according to the 
invention then gives an iterative estimate of the 
discrete channel starting from these initial conditions 

{3 ,(0) }m> * According to the SAGE extension of the EM 
algorithm, the optimum channel estimator according to 
the MAP criterion iteratively calculates the re- 
estimate {3 !(d+l) ]u starting from the estimate {g^}^ 1 , for 
each path in turn. During each iteration, the channel 
estimator calculates the Q 1<d+1) vectors, where 1 = 0, 
1,..., L-l, by successive approximations starting from 

the {g^" 0 }^ and {G m(d) t' vectors calculated beforehand. 
According to the invention, the Jc th component of the 
(d+2) th re-estimate G ltd+1) of the G 1 vector associated with 
the I th path is obtained using the expression: 

Gr=w ] k x 
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where 

w k = j + jsj /^ lk = ^ n -< L - , )'^ n -< L - 2 >'*" ,J ^ n * 

and C™ (d) is the n th component of the reconstruction C m(d) , 
using the KARHUNEN - LOEVE decomposition of the jn th path 
of the discrete channel in the d th iteration. This 
reconstruction is given explicitly by 

N-l 

£m(d) ^^Qm(d)gmk 

k=0 

As before, the weightings w ! k take account of the 
quality of the contribution of the vector B lk to the 
appropriate representation of the discrete channel. 
The higher the contribution of an eigenvector compared 
with the noise variance, the closer its weighting 
becomes to 1 and it plays an increasingly important 
role in the reconstruction of the channel estimate. 

In order to estimate the appropriate 
representation G 1(d+1> of the 2 th path in accordance with 
the previous expression, the receiver starts by 
calculating an estimate 

I-l L-l 

Rv^^mcd+DA Vr m(d) A 

m=0 m=l+l 

of C' n A n _, by subtracting the contributions of the other 
paths m * 1 . This subtraction is made assuming that 

the vector A n = (A n _ (L _ 1} , A n . (L . 2) A n ) is transmitted, 

which is why it is useful to weight this expression 
later by the conditional a posteriori probability 
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P^A n |R 0 {G m(d+1) }^ J of actually having transmitted 

this vector. Once this estimate of C 1 n A n _ 1 has been 
made, the receiver multiplies it by the conjugate A*_, 
of A n _, to become independent of the modulation. This 
operation is done always assuming that the vector A n has 
been transmitted, which is why it is necessary to 
weight the resulting estimate: 

f 1-1 L-l \ 

p Vr m < d+,) A _Vr m(d) A a* 

K n ~~ 7,^n A n-m Z^n ^n-m A n-1 
^ m=0 m=l-l j 

of C[ by the conditional a posteriori probability of 
this vector having been transmitted. Obviously, the 
final estimate 

r. -£cr +,) A n _ m - Zc: (d,A n = m k-. p ^ 

m=0 m=l+l J V J 

of C\ is obtained taking account of the contribution of 
all possible values of A n . 

Once this estimate C ! n of the n th sample of the I th 
path has been obtained for n between I and N+l-1, all 
that remains is for the receiver to project it on each 
of the vectors B lk of the orthonormal base associated 
with the appropriate representation of this path. The 
projection on the Jc th vector B lk is then weighted by 
to take account of the quality of its contribution to 
the estimate of the 1 th path of the channel in the 
(d+l) ch iteration. 

In practice, only a few basic vectors make an 
important contribution to the representation of each 
path. In this case, the number of vectors used in the 
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KARHUNEN-LOEVE transformation and the inverse 
transformation may be reduced significantly with no 
perceptible drop in performances, while guaranteeing a 
significant drop in the complexity. 

An iterative estimate of the discrete channel may 
be made a finite number of times D such that the global 

estimate fc I(D) £! (obtained from {3 I(D) ]£ using the 
KARHUNEN-LOEVE transformation) guarantees that any 
deterioration in the receiver ' s performances compared 
with the optimum solution that would be obtained after 
an infinite number of iterations would be 
imperceptible. This iterative estimate assumes that 
conditional a posteriori probabilities 



are known for all possible 



values of the An = (A n . (L . 1J# A n _ {L . 2) , . . . , A n ) vectors, for 
n = 0, 1,..., N+L-2 . These conditional probabilities 
may be calculated precisely using the BAHL algorithm, 
also known as the BCJR algorithm [6]. They may also be 
calculated approximately and easily using linearized 
variants of the BAHL algorithm [10], 

For example, the following describes an outline of 
the calculation of the conditional probabilities 

P^A n R 0 {G m(d+,) }m^ using the BAHL algorithm. Note 

that according to figure 8, the sample R n obtained at 
instant n depends on the symbol A n transmitted at this 
instant and the L-l symbols A n _ x , A n _ 2 , A n - ( L-i> 
preceding it. Thus, the modulator - multi-path channel 
- demodulator assembly may be modeled using a MARKOV 
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discrete time sequence. This MARKOV sequence is 

described by a mesh, in which the states represent the 
contents of the offset register in figure 8 at each 
instant n and in which the labels of the branches 
represent the non-noisy version of possible samples at 
the output from the whitener filter. For example, 
figure 11 shows a mesh view of a discrete channel 
composed of L=2 paths for a modulation MDP2 with two 
states {-I, +1} and a length of the timeslot N=4 . 

This type of mesh represents the expected output 
from the channel as a function of time and depending on 
the transmitted symbols. In the case shown in which it 
is assumed that the channel has two paths, the channel 
output depends on the symbol transmitted at the instant 
considered and at the previous instant. For example, 
starting from instant n=l, it is possible to transmit 
either -1 or +1, knowing that nothing was transmitted 
at the previous instant 0 (initial node at the left 
end) . If -1 was transmitted, then it would be expected 
that (-lX^o + (°) c o' namely -C° 0 , would be received at the 
output from the channel, which is the quantity that 
appears on the first up branch towards the -1 node. If 
+1 was transmitted, it would be expected to receive 
(+l)Co +OCo , namely C° Q , which is the quantity that 
appears on the first down branch towards the +1 node. 

For n=l, if -1 was transmitted and if -1 had been 
transmitted at the previous instant n=0, then it would 
be expected that (-1) C? + (-1) C] , namely -Cf-Cj would be 
received. If +1 was transmitted, it would be expected 
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that (+1)C? -Cj = C° -Cj would be received. Similarly, if 
+ 1 had been transmitted at the previous instant, it 
would be expected that -C°+C[ would be received if -1 
had been transmitted, or that C^+C} would be received 
5 if +1 had been transmitted. These four values appear 
on the four branches at n=l. The timeslot is completed 
by sending 0, and the expected quantity is -C\ or C\ 
depending on whether -1 or +1 was transmitted at the 
previous instant n=3 . Therefore the mesh actually 
10 represents the quantities expected at the channel 
output at each instant, as a function of the symbols 
transmitted. 

In accordance with this representation of the 
discrete channel, each branch of the mesh determines a 

15 unique value of the vector A n = {A n _ (L _ 1} , A n _ ( ^ 2) , AJ , 

since firstly the first L-l components are given by the 
initial state S n = (A n _ (L _ lt , A n _ (L _ 2} , . .., A n _ x ) of the 
branch, and secondly the last L-l components are given 
by the final state S n+ =(A n _ (L _ 2) , A n _ (L _ 3) A n ) . The 

20 problem of calculating conditional a posteriori 



reduced to a calculation of the probabilities of the 
branches of the mesh depending on the received samples 
vector R and the intermediate estimate of the 
25 normalized discrete channel given by the vectors c"" d+1) 
for the first 1 paths m = 0, 1, 1-1 and the 

vectors C 1 " 1 ^ for the last L-l paths m = 1, 1 + 1,.. . 
L-l. Consequently, the branches of the mesh must be 
labeled by the coefficients of this intermediate 



probabilities 




is 



therefore 
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estimate and not by the coefficients of the real 
discrete channel . 

The BAHL algorithm can be used to efficiently 
calculate conditional a posteriori probabilities of the 
5 branches of the mesh by factoring the expression of the 
probability of the received vector R, depending on each 
of the sequences A that may be transmitted, in order to 
avoid duplicating all operations common to transmitted 
sequences with the same past or the same future at a 
10 given instant. Consequently, this algorithm uses the 
a n (S) and P n (S) functions that give a preview of the 
past and the future respectively for each instant 
n=0, 1, N+L-l, and each possible value S of the 

state S n . 

15 The function Ot n (S) is the joint probability of 

being in state S n =S at instant n and having received the 
samples {R j . This function is calculated by 

factorization of the calculations for sequences with a 
common past using the recurrence: 

20 a n (S) = Xa„-,(S')Y„(S\S) 

S" 

where Y n (S',S) is the function that gives the 
probability, depending on an initial state of the mesh 
S^^S', of changing from state S n =S and receiving the 
sample at the same time. This expression shows that 

25 information about the past may be acquired either 
through probabilities associated with the instant 
immediately preceding the current instant or through 
the probabilities of transitions of the mesh branches. 
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The significance of the P n (S) function is similar 
to the significance of the oc n (S) function. It is the 
probability of having received the samples {R-klkln" 2 
depending on the initial state S n =S at instant n. This 
function may also be calculated by successive 
approximations starting from the future by the 
recurrence 

P„(S) = XP„ +1 (S)Y„ +1 l(S\S)J 

S' 

Here again, this recurrence proves that knowledge 
of the future for a state S n =S may be obtained starting 
from the characteristics of the states S n+1 =S ■ 
immediately following this state and the branches of 
the mesh joining them together. For a branch of the 
mesh characterized by the initial state S n = S' = (A a _ a _ 
i>' A n-u- 2 >' A n-i) and the final state S n _=S= (A n _^ 

x ) , . . . , A n-1/ A n ) , the Y k (S',S) function is given by: 



Yn+1 (S\S) = P(A n )-i-exp 



f 



1 

N 0 



1-1 L-l 



R.-Zcr->A M .-xcr d) A M 



m=0 



2 >t 



where P(A n ) is the a priori probability of sending the 
symbol A n . These a priori probabilities may 

characterize a possible coding of the actually 
transmitted modulated symbols. They may also 

characterize the reference symbols that are known to 
the receiver and that have only one possible value. 
For these reference symbols, the a priori probability 
P(A n ) is equal to 1 for the value actually transmitted 
and 0 for other possible values of the modulation. 
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The following initial values have to be used so 
that recurrences can be made on the a n (S) and P n (S) 
functions: 

a 0 (S0 = {0, 0 0)) = 1 and P a . x _ l = (S^ = (0, 0 0)) = 1 

Starting from a calculation of the ot n (S) and P n (S) 
functions using the BAHL algorithm, it is easy to 
obtain the probabilities for branches of the mesh, 
depending on the received samples and the intermediate 
estimate of the discrete channel using the formula: 

oc„(S , )Y n+1 (S',S)P n+I (S) 



p(a d =A|R,{G m <^ 



0 O (S O =(O,O,...,O)) 

where A n = (A a _ (L _ 1} , A n _ {Zj _ 2} , . AJ , S' = (A n _ (L _ 1} , A a , (L , 2l , 

. . . , A n _ 2 ) and S = (A n _ (L _ 1} , . . . , A n _ lf AJ . 

However, note that the BAHL algorithm is fairly 
complex and that in practice it is better to use 
approximate algorithms that can significantly reduce 
the complexity of the calculations without excessively 
reducing performances. The use of these approximate 
algorithms is recommended particularly when the number 
of paths and/or modulation states is high. 

After D iterations, the receiver has an estimate 

of the direct channel jc l(D) ]|lo close to optimum value 
depending on the maximum a posteriori criterion that is 
normally obtained after an infinite number of 
iterations. Starting from this channel estimate and 
samples received during the timeslot to be processed, 
the receiver may output either firm decisions using the 
VITERBI algorithm, or weighted decisions using the SOVA 
algorithm or the BAHL algorithm. 



SP 16348. C RS 



47 



The VITERBI algorithm searches for the most 
probable transmitted sequence. Therefore, it can 

reduce the probability of decision errors on these 
transmitted sequences. The SOVA algorithm is a VITERBI 
algorithm with weighted outputs. It enables a more 
precise characterization of the reliability of 
decisions made with the VITERBI algorithm. The BAHL 
algorithm can determine the a posteriori probability 
density for each transmitted symbol, depending on the 
received samples and the obtained channel estimate. 
Therefore, it can determine the most probable value of 
each transmitted symbol and therefore reduce the error 
rate. Weighted outputs are obtained through symbol 
probabilities depending on received samples. 

In all cases, weighted outputs are only useful if 
corrector decoding or source coding with weighed inputs 
is used in order to increase the receiver performances. 

For example, within the framework of the BAHL 

algorithm, the receiver uses the {c^}^ 1 estimate of the 
discrete channel to calculate the ot n (S) and (3 n (S) 
functions. These functions can then be used to 

determine the probability distribution for each symbol 
A n using the formula: 



A = A 





f 




A n =A 


j A 





R.{3 ,CD) l2 = X P A n=AR.{G ,(D t 



J 



All the operations for iterative estimation of the 
discrete channel and detection operations described 
above are shown in figure 12 . The channel estimator 

120 uses means defining the {b'^J bases, 1=0, 1, 
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L-l, and means 122 calculating the weighting 
coefficients {w ,! k E and {w^}^ . The estimator 120 uses 
a first coarse estimate of the channel starting from 
the elementary reference sequences alone by circular 
correlation with the corresponding CAZAK sequences. 
This first estimate is symbolized by the block 3 0 
reference 12 01 that outputs the appropriate 

representation of the initial estimate {G'^jjlo . The 
estimator 120 then uses an order d re-estimate, 
represented by block 3 d , reference 12 02 that outputs 
the appropriate representation of the channel estimate 

fc I(D) jilo ' followed by an order d+1 re-estimate 
represented by block 3^ reference 12 03 that outputs 
the representation of the channel estimate {g 1 ^" 0 }^ . 
Finally, the estimator 120 uses a final order D re- 
estimate represented by block 3 D reference 12 04 that 
outputs the appropriate representation of the channel 

estimate {g I(D) IJIq 1 . The estimator uses this appropriate 
representation of the estimated channel to output the 
final estimate {c ,(D) }Ilo of the discrete channel using 
means 1205, and using the KARHUNEN-LOEVE 

transformation . 

Firm or soft decisions on the transmitted symbols 
are supplied to an error corrector decoder or a source 
decoder by means 123 that may use the VITERBI algorithm 
[3], the SOVA algorithm [2], the BAHL algorithm [6] or the 
linearized BAHL algorithm [10]. 
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According to the invention, the iterative re- 
estimate of the channel is made D+l times and is used 

to calculate {g 1 ^ , {G ,(1) fo {G^to successively. 

Figure 13 illustrates the L steps making up each 
iteration and that are necessary for successive 
estimate of each of the L paths. The (d+l) th iteration 
3 J , of the channel estimator uses the means 121 and 122 

d+l 

shown in figure 12 (providing bases {b"*}^ and 

weighting factors {w™ jj^J ) , carries out a first step 
based on the BAHL algorithm and represented by the BAHL 0 
block reference 13 0, to produce conditional a 



posteriori probabilities P 



R.{G m(d, tl 0 ' 



starting from 



n i- -»m— u 

V J 
the cf h re-estimate of the channel. Starting from these 
conditional probabilities, the (d+l) th iteration then 
executes a second step, based on the SAGE algorithm and 
represented by the SAGE 0 block reference 131 to enable 
the (d+l) th re-estimate G 0<d=1) of the first path. 

Starting from the (d+l) ch re-estimate of the 1-1 
first paths, the (d+l) th iteration executes a third 
step, once again based on the BAHL algorithm and 
represented by the BAHL block, reference 132 to produce 
conditional a posteriori probabilities 

P^A n |R.{G m(d+I) £;o. ,{G m(d) E J starting from the d* h re- 
estimate of the L-l last paths of the channel and the 
(d+l) th re-estimate of the first 1 paths of this 
channel. Subsequently, the (d+l) th iteration executes a 
fourth step based on the SAGE algorithm and shown by 
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the SAGE 1 block reference 133, to enable the (d+1) re- 
estimate G 1(d+1) of the I th path. 

Finally, starting from the (d+l) th re-estimate of 
the first L-l paths, the (d+l) th iteration executes a 
5 final step based on the BAHL algorithm and represented 
by the BAHL 1 _ 1 block reference 134, to output conditional 



starting from the <? h re-estimate of the last channel 
path and the (d+l) th re-estimate of the first L-l paths 

10 of this channel. The (d+l) th iteration then executes a 
final step, based on the SAGE algorithm and represented 
by the SAGE L _ 1 block reference 13 5, to make the (d+l) ch 
re-estimate G (tj " 1> ' <d " 1) of the L th path. 

The first iteration symbolized by block 3 0 

15 reference 1201 in figure 12, cannot use the BAHL and 
SAGE algorithms to calculate the first estimate of the 
discrete channel. According to the invention, the 
reference symbols for which a priori and a posteriori 
probabilities are well known to the receiver are used 

20 to make this first estimate. This first estimate is 
made by correlating the received samples depending 
solely on these reference symbols, with the 
corresponding CAZAC sequences. It uses means 122 in 

r ,| tn-1 

figure 12 to calculate the weightings tw k / kai0 . 




a 



posteriori 



probabilities 
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